Efficiency of dynamical decoupling sequences in presence of pulse errors 
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For a generic dynamical decoupling sequence employing a single-axis control, we study its efficiency 
in the presence of small errors in the direction of the controlling pulses. In the case that the 
corresponding ideal dynamical-decoupling sequence produces sufficiently good results, the impact 
of the errors is found to scale as with negligible first-order effect, where ^ is the dispersion 
of the random errors. This analytical prediction is numerically tested in a model, in which the 
environment is modeled by one qubit coupled to a quantum kicked rotator in chaotic motion. In 
this model, with periodic pulses applied to the qubit in the environment, it is found numerically 
that periodic bang-bang control may outperform Uhrig dynamical decoupling. 
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I. INTRODUCTION 

Dynamical decoupling (DD) has attracted lots of at- 
tention in the past years, due to its potential applica- 
tion in quantum information processes. The basic idea 
is to insert a sequence of controlling pulses within a time 
period of interest, such that the system of interest can 
be effectively decoupled from its environment. As a re- 
sult, environment-induced decoherence can be effectively 
suppressed and initial coherence in the system can be 
preserved. 

Several DD schemes have been proposed. For example, 
the so-called periodic "bang-bang " control [ll-Q can sup- 
press decoherence up to the order of 0{T^) for a given pe- 
riod T of coherence preservation, while the Carr-Purcell- 
Meiboom-Gill sequence has an 0{T^) efficiency 0, i|. 
A better result can be obtained by an approach called 
"concatenated dynamical decoupling" d, Q , with an ef- 
ficiency of the order of ©(T^"*"^) achieved by 2^ pulses. 
Recently, a remarkable progress has been made by Uhrig, 
showing that decoherence of a single spin, which is in- 
duced by a spin-boson bath, can be suppressed up to 
the order of 0(T^+^) with only N pulses Q. Later, the 
Uhrig dynamical decoupling (UDD) was conjectured 
and rigorously proved to be model-independent for the 
pure dephasing case ^ . More recently, it has been found 
that different types of UDD sequences may be integrated 
to obtain better results [lOj and, with appropriate exten- 
sion, to work for a system with two spins as well [Tl| . 
Meanwhile, the efficiencies of UDD and its generaliza- 
tion has been beautifully demonstrated experimentally 

[TMl 

In practical application of a DD sequence, an impor- 
tant topic is its robustness, i.e., the influence of non- 
idealness of the controlling pulses in the efficiency of the 
DD. The non-idealness may come from finite widths of 
the pulses and/or small deviation of the actual direc- 
tions of the pulses from their designed directions. In the 
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case that the non-ideal parts of the pulses with finite 
width satisfy certain symmetry requirements, a general- 
ized UDD can be found [9|]. However, in a generic cir- 
cumstance, accumulation effect of small imperfections in 
the pulses may have significant consequences when the 
number of pulses is not small. There have been several 
investigations in the impact of systematic pulse errors for 
some specific DD sequences, both experimentally [l5l - [l8j 
and theoretically [17-2Q|. 

In this work, for a generic DD employing a single-axis 
control, we study accumulation effects of small random 
errors in the direction of the controlling pulses. Analyti- 
cal derivations are given in Sec|TTl In particular, we show 
that a DD, which has a sufficiently good performance in 
the ideal case, is robust in the sense that the pulse errors 
have negligible first-order effect. 

In numerical investigation, we employ a model in which 
there are periodic pluses in the environment, which are 
not directly applied to the qubit of interest. Using this 
model, we can test our analytical predictions, and also 
study another topic of interest, namely, the influence of 
high-frequency cutoff in the efficiency of UDD. When the 
spectrum of the environment has a sharp high-frequency 
cutoff, UDD has been found outperforming all other 
known DD sequences and is regarded as optimal. How- 
ever, for environments with soft cutoffs in the spectra, 
there is no reason to expect that UDD is optimal; in 
fact, recently it has been shown that protocols with pe- 
riodic structure, such as the Carr-Purcell-Meiboom-Gill 
sequence, ma y h ave a better performance for this type of 
environment j2Tl - [2l | . 

Specifically, In Sec lIIIl we study a model, in which the 
environment is simulated by a second qubit coupled to 
a quantum kicked rotator in chaotic motion. Previous 
study shows that this model, though a single-particle dy- 
namical system, may simulate a pure-dephasing many- 
body bath (Caldeira-Leggett model [11]), as well as some 
non-Markovian environments (26| . Our numerical simu- 
lations show that for this type of environment the peri- 
odic bang-bang control may outperform UDD. Finally, 
conclusions and discussions are given in Sec lIVI 



2 



II. ANALYTICAL STUDY OF IMPACTS OF 
SMALL PULSE ERRORS 

In this section, we first recall essential properties of 
DD, then, derive expressions for the influence of small 
pulse errors in the efficiency of a DD that employs a 
single-axis control. 



A. Dynamical decoupling scheme 

Let us first recall some essential properties of a generic 
DD. Consider a qubit S and its environment £, under a 
Hamiltonian 

H = Hs + Hi + He, (1) 

where Hs and Hg are the Hamiltonians of S and £, re- 
spectively, and Hj is the interaction Hamiltonian. The 
self-Hamiltonian Hs is assumed to be a constant, which 
can be taken as zero, Hg — 0. We further assume that 
the interaction Hamiltonian Hj is commutable with cr^, 
[az,Hi] = 0, where cr^ is the ^-component of the Pauli 
operator for the qubit S. 

Suppose the qubit S lies initially in a state, which 
is not an eigenstate of the interaction Hamiltonian Hi. 
Schrodinger evolution under the Hamiltonian H, given 
by a unitary operator U{T, 0) for a time period T, may 
induce decoherence to the state of the qubit S. The pur- 
pose of a DD is to preserve coherence in the initial state 
of S within the time period T by inserting a sequence of 
pulses, e.g., TT-pulses for the x direction. Using ti = SiT 
of i = 1, 2, • • • , n to indicate the instants at which pulses 
are applied, with totally n pulses, the time-evolution op- 
erator is now written as 

R=U{T,tn)a.,Uit„,t„.i)a,---Uit2,h)a,U{h,0). (2) 

To measure preservation of coherence in an initial state 
in the x-direction, one may consider measurement on the 
observable ax , which gives the signal 

s{T) = (t \Dl{-n/2)R^<JxRDy{~TT/2)\ t), (3) 

where | f) indicates an eigenstate of az and Dy{—TT/2) 
rotates it to an eigenstate of a^- After some derivation, 
one gets 

s{T)^mc{i\R^axR\t)- (4) 

The coherence is perfectly preserved if s(T) — 1. 

In UDD, TT-pulses are applied to the qubit S at times 
tj = (5jr, where 

5j = sin2[7rj7(2n + 2)] (j = 0, 1, 2, • • -n, n + 1). (5) 

For UDD, s(r) = 1-0(T2"+2). In a periodic bang-bang 
control of DD, pulses are applied at times with 8j = j /n. 



B. s(T) expanded to the second-order term of error 

As discussed in Introduction, controlling pulses in a 
DD may be subject to small random errors in its direc- 
tion. Let us consider small random deviation in the y 
direction for 7r-pulses in the x direction. In this case, ax 
in the ideal time-evolution operator in Eq. ^ should be 
replaced by 

CTe = Exax + SyCTy, (6) 

where e^, = y^l — and Ey is a small random number 
with Gaussian distribution. 

Here, ^ is the dispersion of the random number, with 
^ ^ 1. Then, the time evolution operator, denoted by 
_Re, is written as 

Re = U{T,tn){^JT~^el^ax + Ey^nf^y) 

• • • U{t2M){^l~el^ax + ey^iay)U{h,Q), (8) 

where we use Sy^i to indicate the value of £y for the i- 
th pulse. Now, the signal, denoted by Se{T), has the 
following expression, 

Se{T)^mc{i\RlaxR,\t)- (9) 

To evaluate Se^T), we substitute Eq. ^ into Eq. dH) 
and expand the result in the power of Sy. Up to the 
second-order term of ^, we write the signal in the follow- 
ing form, 

s,{T) = s,o(r) + sei(r) + s,2{T) + O(C'). (10) 

Definitions of the first three terms on the right hand side 
of Eq. ([T0|) will be given below, when they are treated 
separately. 

The first term on the right hand side of Eq. (fTO|) is 
obtained by considering only the contribution of Sxax in 
each de, 

s,o{T) = mc{i\RloaxReo\t), (H) 

where 

ReO = U{T, tn)^l - el„axU{tn, tn-l) 

n 

=n(i-<^)^^- (12) 
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Here, R is the time evolution operator for the case of ideal 
pulses in Eq. ©. Substituting Eq. dH]) into Eq. (HT]), we 
find 



k=l 



It is seen that, for sufficiently small ^, 

s,o{T)-s{T)^~nes{T). 



(13) 



(14) 



Therefore, Sei(r) gives a first order correction (~ 
to the ideal signal. 

The first order correction Sei(r) also depends on the 
quantity q{T). To give an estimation to q{T), we note 
that s{T) +iq{T) = {i \R^C7^R\ t), [see Eqs. Q and 
([22| ]. Hence, 



Thus, deviation of Sgo{T) from the ideal signal s{T) is of 
the order of {^^/n)'^. 

Next, we calculate the second term on the right hand 
side of Eq. pUj) . which is the contribution of those multi- 
plication terms that include only one (SyCy) term in each 
of them. Noticing that, up to the second-order contribu- 
tion of Sy, Ex in this second term can be taken as 1, we 
have the following expression for it, 

Sei(T) = fHe(i \Rl^axReo + Rlo<yxRei\ t), (15) 
where 

= U{T,t U ■■■U {t2,ti)<JxU{ti,0) 

+ U{T, tn)(TxU{tn,'tn^l)£y,n-l(^y ' ' ' U {t2,ti)axU ih^^) 

H 

+ C/(T, tn)<TxU{tn, tn-l)(Jx ■ ■ ■ U {t2, h)ey sCTyU {h, 0) . (16) 

Making use of the fact that ay — iax<Jz and [gz, H] = 0, 
we find 



\q{T)\ = ^\{i\R^axR\^W-s^{T) < ^1^^, 

(23) 

where we have used the fact that |(| {R'^axRl t)l ^ 1- 

Finally, we discuss the third term on the right hand 
side of Eq. (fTO|) . which includes all multiplication terms 
that have only two {SyCTy) terms, 

Se2{T) = 5nc(; \Rl^axRel + RU<yxRe2 + Rl2<JxReo\ t) , (24) 

where 

Re2 — U{T,tn)ey^n<^yU{tn,tn_i)ey^n-l<^y 
)'^x ■ ■ ■ U{t2,ti)axU{ti,0) 
+ U{T,tn)Sy^nO'yU{tn,tn-l)<Tx 

U(tn,tn-l)ey^n-2<yy ' ' ■ U {t2,ti)axU {ti,Q) 

+ ^U{T,t„)axU{tn,tn-i)ax 

■ ■ ■ U (^3, t2)£y,2CryU{t2, tl)eyAO'y (7(ti,0). (25) 

Following a procedure similar to that for -R^i, we find 

1 



k=lj{<k) 



e2^ 2^ 2^ (-l)J+^-+^ej,,,£,,fci?,0 - ±-r'^^e0.(26) 



Rel = lY.^-lf^^£y,kReQ(Tz. 



(17) 



k=l 



Substituting Eq. ^ into Eq. (d?]), then into Eq. ([T5|). 
we obtain 



2^^(-l)^+ley^fc(i \Rl^axRe^\ t) 



fc=i 



,(18) 



Typically, one has the following estimate. 



k=l 



hence. 



,l(T) ^ ±2Cy^^Z \RloCTxReo\ t) 



(19) 



(20) 



Substituting the above obtained expressions of R^o, 
and Re2 into the expression of 562(7") in Eq. ((24l) . 
after some derivation, we obtain an expression for 5^2 (T) . 
Then, making use of results obtained above in Egs. p^ 
and (|21l) . finally, we find 

s,iT) - s(r) ^ ±2q{T)CV^ + C2en + O(e'), (27) 

where IC2I is of the order of 1. 

Of particular interest is the case for a DD with good 
performance, i.e., with s(T) ^1. In this case, the in- 
equality (j23p shows that q{t) is small. In particular, in 
the case that \q{t)\ <^ ^y/n, the first order term on the 
right hand side of Eq. ([27|) can be neglected and we have 



(28) 



,{T)~s{T)^C2en + 0{e), 



Making use of Eq. p2|). after simple algebra, it is found 
that 



where 



S6i(T)~±2eV^g(r), 
q{T) = Jm{i\R^axR\ t)- 



(21) 
(22) 



scaling as ^^n. For example, for a UDD with s(T) — 
1 - 0(T2"+2), q{T) is of the order of 0(T2"+2) or less, 
hence, for a fixed ^, |?(r)| ^ for a sufficiently large 

In concluding this section, we remark that it is 
straightforward to generalize the above discussions to the 
case of pulses with more generic random errors in their 
direction. In fact, in this generic case, the time-evolution 
operator can be obtained by replacing ax in Eq. ([2]) by 



(29) 
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0.0 0.5 1.0 0-0 0-5 1-0 

t/T t/T 

FIG. 1: (Color online) Variation of the signal s{t) in Eq. p4|) 
with time for ideal UDD and ideal periodic bang-bang control. 
Both UDD and periodic bang-bang control have better perfor- 
mance for larger values of n, the number of controlling pulses. 
Results of UDD (circles) are better than those of the periodic 
bang-bang control (solid curve) for n = 50. But, for larger 
values of n, the periodic bang-bang control outperforms UDD. 
Parameters are lja = 1.5 x 10^, A = 10^, g = 100, = 2^^ and 
T — 50To, and the kicked rotator is in the chaotic region with 
K = 10^ The same parameters are also used in the following 
figures. 



where Sx = y 1 — — and both Sy and Sz are small 

random numbers with Gaussian distribution. We have 
found results that are qualitatively the same as those 
discussed above for the case of in Eq. In particu- 
lar, we have found similar estimations as those given in 
the relations ^ and (|28l). 



III. NUMERICAL SIMULATIONS 

In this section, we discuss numerical simulations we 
have performed for two purposes. One is to check ana- 
lytical predictions given in the previous section, the other 
is to study the influence of kicks in the environment in the 
efficiency of UDD. In fact, since instant kicks in the envi- 
ronment may have non-negligible high-frequency compo- 
nents, the performance of UDD for such an environment 
may be not so good as that for an environment with a 
sharp high-frequency cutoff. 



A. The model 

We consider a model, in which there is a qubit S of 
interest and an environment £ that is composed of a sec- 
ond qubit A and a quantum kicked rotator denoted by 
B. The qubit S has interaction with A only, while A 
interacts with both S and the kicked rotator B. The 
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FIG. 2: (Color online) Dependence of l — Se{T) (solid squares) 
on 5^ for a UDD sequence, with n = 500 and ^y/n <^ 1. 
The solid curve is a straight line, indicating a linear depen- 
dence, in agreement with the prediction of Eq. H28p . For 
the corresponding ideal UDD, 1 — s(T) ~ 10~*. Since 
1 — Se(r) 2> 1 — s(r), the infiuence of the pulse errors is 
profound. 

Hamiltonian is written as 

H = Hs + Ha + Hb + Hsa + Hab, (30) 
where the self-Hamiltonians are Hs — 0, Ha—uja<^x^ and 

ffs = y + fccos0^(5(i-jTo), (31) 

i 

with To the period of kicking. The interaction Hamilto- 
nians are Hsa = 5f f ® and 

Hab = Xa^ cos eJ2Kt- jTo). (32) 

3 

Here, for clearness we write explicitly the superscript S 
in the Pauli operator for the qubit S. 

This model has been studied in Ref. [26], showing 
that the kicked rotator, though a single-particle dynam- 
ical system, may simulate a pure-dephasing many-body 
bath (Caldeira-Leggett model [25|), as well as some non- 
Markovian environments. Here, we are interested in the 
chaotic region of the kicked rotator, to simulate some 
random properties of the environment. 

The time evolution for one period Tq is given by the 
unitary operator 

Uj,^ = e-»('^ATorT^-Hro3rTf 0f7^)g-iToJ^g-j(fe+A£7^)cose^(-33^ 

where H has been set unit. An effective Planck constant 
can be introduced, ftoff — Tq = 2tt/N, where N is the di- 
mension of the Hilbert space of the kicked rotator. The 
classical limit is obtained by letting Tq — >■ and A: — >■ oo 
while keeping K = kTg fixed, and the classical counter- 
part is defined on a torus [0,27r) ig) [0,27r). The kicked 
rotator has a chaotic motion for K larger than 6 or so. 
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FIG. 4: Variation of the signal s^{t) with time for a UDD 
with non-ideal pulses, ^ = 10~^. The UDD with n = 500 
(right panel) gives results worse than those of the UDD with 
n = 200 (left panel). 



FIG. 3: Variation of 1 — Se(r) with C\/n for ^ UDD sequence 
with n = 500. 

B. Numerical results 

Let us first discuss the performance of ideal UDD and 
ideal periodic bang-bang control in this model. We are 
interested in the case that T ^ Tq. In this case, due to 
the time-dependency of the Hamiltonian of the environ- 
ment, Uhrig's strategy for deriving 1 — s{T) ~ 0(r-^"+^) 
is not applicable. It is of interest to know whether UDD 
is optimal or not in this case. 

In numerical simulation, we calculated s{t) defined by 

s{t)^D\z{i\R\t)a,R{t)\t), (34) 

where R(t) is obtained by truncation of the time evolu- 
tion operator R in Eq. ^ at an intermediate time t < T. 
Figure [T] shows that both UDD and periodic bang-bang 
control have better performance with increasing number 
n of the controlling pulses. For n = 50, UDD outper- 
forms the periodic bang-bang control in the whole time 
region (0,T]. However, with increasing n, results of the 
periodic bang-bang control become better than those of 
UDD for t > T/2. 

Next, we check analytical predictions given above for 
the impact of random errors in the direction of the con- 
trolling pulses, in particular, the behavior of Se{T) in 
Eq. (|28l) . We use the general form of in Eq. (f29| , with 
the same dispersion ^ for Sy and . An example is given 
in Fig 121 showing linear dependence of 1 — S£(T) on To 
check details of agreement with analytical predictions, we 
have numerically computed the corresponding ideal UDD 
and found that it has a good performance with 1 — s(T) ~ 
10"''. This gives \q{T)\ < 10" ^ [see the estimate in the 
inequality Hence, \q{T)\ < ~ 10" ^ for solid 

squares shown in the figure, as a result, Eq. (1^51) should 
hold with s(T) - s^iT) ~ 1 - Se(T). Thus, results in 
Figl2] indeed confirm the prediction of Eq. (pS)) . 

We have also numerically studied the case of large ^\/n- 
In this case, higher order terms of ^ Eq. (I?71) should be 
considered. Anyway, this expression of s^{T) suggests 



that 1 — Se{T) may be large for f^n ^ 1. Indeed, numer- 
ical results support this expectation, as shown in Fig. |31 
Figure [2] shows that for a fixed n, the influence of pulse 
errors may become large with increasing ^, such that 1 — 
Se(T) 3> 1 — s(T). In Fig|31 we show that for a fixed error 
dispersion ^, deviation of s^it) from 1 enlarges when n 
is increased, where Se(i) is defined by Eq. (p4| with R{t) 
replaced by the corresponding Re{t). 



IV. CONCLUSIONS AND DISCUSSIONS 

In this paper, we have analytically studied the effi- 
ciency of a generic DD with a single-axis control, when 
the controlling pulses are subject to small random errors 
in their direction. We have derived expressions for the 
influence of the pulse errors up to the second-order term. 
When the ideal DD has a sufficiently good performance, 
the influence has a negligible first-order effect; in this 
sense, good DD are relatively robust. 

We have tested the above analytical predictions nu- 
merically and shown that accumulation of small pulse 
errors may have significant influence in the efficiency of 
DD. For an environment with kicks applied on some part 
of the environment, it has been found that the periodic 
bang-bang control may outperform UDD. 

A natural question would concern the possibility of 
having negligible first-order effect of pulse errors in more 
generic situations, e.g., in the case of more than one lay- 
ers of controlling pulses with different single-axis control 
in different layers. For pulse errors appearing only at 
certain fixed layer with single-axis control, results of this 
paper may be generalizable. However, the more generic 
situation with pulse errors in different layers, as well as 
the case with finite width of the pulses, are much more 
complex and further investigation is needed before a def- 
inite conclusion can be drawn. 
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